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Abstract. We prove some regularity results for a class of two dimensional 
non-Newtonian fluids. By applying results from [Dashti and Robinson, Non- 
linearity, 22 (2009), 735-746] we can then show uniqueness of particle trajec- 
tories. 

1. Introduction 

In this paper we consider the following system of partial differential equations 

(1.1a) u t -i^)Ait-i/idiViS , (7>u) + (u-V)ti + V7r = / in [0,T]xO, 

(1.1b) divu = in [0,T] X CI, 

(1.1c) u(0) = u in Q, 

where f2 denotes either a two-dimensional bounded domain or the two dimensional 
flat torus, the vector field u — (u\, U2) is the velocity, the scalar n is the kinematic 
pressure, the vector / = (fi, f%) is the external body force, uq is the initial velocity, 
and i/o, v\ are positive constants. We denote by 

Vu := i(Vu + V/) = ^(djUi + dittj) for i,j = 1, 2, 

the symmetric part of Vtt, the convective term is (u ■ V) U := Y^k=i u kdkU, and S 
denotes the extra stress tensor, defined by 

(1.2) S{Vu) := {5+ \Vu\) p ~ 2 Vu, p€[l,2), 

where S is a non- negative constant. System (11.11) describes a shear-thinning ho- 
mogeneous fluid and for an introduction to the mathematical theory see Malek, 
Rajagopal, and Ruzicka We mainly study the problem, endowed with homo- 
geneous Dirichlet boundary conditions 

(1.3) U| r =° where r = <9ft, 

but we give some remarks also on the periodic case. 

The main goal of this paper is to study the problem of uniqueness for the particle 
trajectories (or characteristics), which are solutions of the following Cauchy problem 

f dX(t) , , , , , , 

(1.4) hr^^^ te[0 < T] ' 

where u is the fluid velocity in (ll.ip . For the 3D Navier-Stokes equations the 
problem of existence of particle trajectories and Lagrangian representation of the 
flow started with the work of Foias, Guillope, and Temam (T3], and related results 
of regularity in W 1 are proved in Chemin and Lerner by means of Littlewod- 
Paley decomposition. The question of uniqueness has been recently addressed by 



Date: March 8, 2013. 



1 



2 



L. C. BERSELLI AND L. BISCONTI 



elementary tools and in a more general context in Robinson et al. [TU1 [2D [22] and 
it is strictly related with uniqueness for linear transport equations. We consider 
here the same problem, in the case of shear-thinning fluids, described by 
To this end, we will study certain regularity properties of the solutions of (jl.ip . 
investigating when the velocity will verify the appropriate hypotheses for uniqueness 
results. 

In particular, classical results concerning Lipschitz continuous fields u (which 
generally can be verified checking that Vu is bounded in the space variables) are 
not easily applicable here, since such a regularity is very difficult to be proved, 
even in the two dimensional case, for ([l.ip . We recall that, restricting to the two 
dimensional case, some C 1,7 -results are obtained in Kaplicky, Malek, and Stara [141 
[T5] in the stationary case. Early results in the time dependent case (but not up-to- 
the-boundary) are those by Seregin [23j . while results in the space-periodic time- 
dependent case have been obtained in [16] . We observe that essentially all the above 
results require the extra-stress tensor S to be slightly smoother than that in (|1.2p . 
In particular, it is requested that the stress-tensor is replaced, for instance, by 
S(Vu) = {S+lVul^^Vu. In any case we study the regularity up-to-the-boundary 
with non-smooth initial data and our results, proved in an elementary way, are 
original. Moreover, the difficulties appearing in the 3D case seem completely out 
of the current mathematical knowledge for such fluids, and this explains why we 
restrict to the two dimensional case. 

Since we want to have elementary proofs (in order to possibly extend the results 
to the widest possible class of solutions and stress-tensors) we will work with the 
classical energy-type methods. Concerning uniqueness of particle trajectories, there 
have been some recent improvements, strictly related with the Osgood criterion 
and with Log-Lipschitz properties of Sobolev functions W s+1 ' q (M. d ) in the case of 
limiting Sobolev exponents such that q = —, In particular we will use the result 
below, proved in [TUJ Theorem 2.1]. 

Theorem 1.1. Let fi be either the whole space WL d , d > 2, a periodic d-dimensional 
domain, or an open bounded subset ofM. d with a sufficiently smooth boundary. Let 
assume that for some p > 1 

u e L P (0,T; W^> 2 {n)) and \ftu € L 2 (0,T; W^' 2 (Q)), 

with u\r — 0, when Q is a domain with boundary. Then, the Cauchy problem (|1.4j) 
has a unique solution in [0, T\. 

The latter result shows that certain (slightly weaker than C 1,7 ) results of Sobolev 
space-regularity can be used to obtain uniqueness for (jl.4[) . On the other hand, the 
W 2,2 (M. 2 ) regularity for fluid with shear-dependent viscosity is another non-trivial 
task (while in 3D proving u € 1^ 5 / 2,2 (R 3 ), seems at the moment out of sight). Some 
recent results (in the stationary case) for second-order space-derivatives appeared 
in [21 IH E] even if the square integrability of second order derivatives is not reached 
in general domains, or if certain limitations on the smallness of the force are not 
satisfied. For the non-stationary case, we recall the result in the space periodic 
setting (obtained uniformly in S > 0) from [S] 111] . 

We also point out that one of the main technical obstructions is represented 
by the pressure and the associated divergence-free constraint. In the case of the p- 
Laplacian systems, in fact, the recent results in Beirao da Veiga and Crispo [3] show 
that u £ W 2 ' q {Q), for arbitrary q, if / is smooth, and under certain restrictions on 
the range of p € (1,2). These latter results are proved in the stationary case, they 
have no counterpart for the p-Stokes system, and most likely they can be adapted 
to the time-dependent case. 
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We point out that in the case of non-Newtonian fluids many features of the 
problem are critical: The type of boundary conditions, the range of p, and if the 
parameter 6 is strictly larger than zero. We will discuss later on some of the 
technical issues of the problem and we will explain why we have to reduce to the 
2D case with i/q,5 > 0. We start by considering the easier case of the periodic 
setting where il is the flat 2D torus T 2 := R 2 /27rZ and we will prove the following 
result. 

Proposition 1.1. Let vq > 0, S > 0, and p G (1,2]. Let be given uq G L 2 (T 2 ) such 
that divuo = and f G L 2 (0, T; L 2 (T 2 )). Then, weak solutions to (jl.l[) satisfy 
\ftu G L 2 (0,T;W 2 ' 2 (T 2 )) and hence problem (|1.4j) admits a unique solution. 

We emphasize that the assumption vq > is crucial in our method. When vq = 
it is possible to prove a regularity result that, although it is not useful to get an 
application of Theorem ll.il seems interesting by itself. See Prop. ETTl cf. Kost [T7] . 

In the Dirichlet case the problem of regularity is more delicate. We will consider 
problem (jl.ip in a domain with flat boundary. We first prove a regularity result, 
by using techniques similar to those used in 8 and formerly introduced, for the 
case p > 2, in [T]. With smooth data, we have the following result. 

Proposition 1.2. Let 6 > 0, v > 0, p G [|,2], u G W 2 > 2 {Vl) n W 1,2 (O) with 
divit o = 0, and f G W rl ' 2 (0, T; L 2 (tt)). Then, weak solutions to Problem ([TTT ]) -<[TT3 |) 
satisfy 



where C depends onp, 5, isq, i>\, ll/lliv^fo, t ; l 2 ); 1 1 i^o 1 1 2,2 7 T. and fl. 

Some hypotheses can be relaxed, since the time regularity is unnecessary for 
the proof of uniqueness of particle trajectories, but the arguments used to prove 
Proposition 11.21 will play a fundamental role to demonstrate our main uniqueness 
criterion for the problem (jl.4[) . The main result of this paper reads as follows: 

Theorem 1.2. Let 6 > 0, v > 0, p G [|,2], u G L 2 (£l) with divw = 0, such that 
(uq ■ n)|p = 0. and f G £ 2 (0, T; L 2 (fl)). Then, weak solutions to (|l.ip - (|1.3p satisfy 
\/tu G L 2 (0, T; W 2 ' 2 (tt)), and consequently (|1.4[) admits a unique solution. 

Plan of the paper. In Section [5] we introduce the notation and we give some 
preliminary results. In Section [3l we consider the space-periodic setting and we 
prove Proposition ll.il Thereafter, in Section 31 we prove a preliminary space-time 
regularity result for the solutions of (ll.ll) - (|1.3p and then we demonstrate Proposi- 
tion 11.21 Finally, in Section [5l we give the proof of Theorem 11.21 



Let us introduce the notation related especially to the problem (ll.ip with Dirich- 
let boundary conditions. The needed assumptions or changes for the space periodic 
case are specified in Section [31 

Throughout the article, when Q is a bounded domain with boundary, it will be 
a two dimensional cube ft =] — 1, 1[ 2 and we denote by T the two opposite sides in 
the X2 direction 

r := {x = (xi,x 2 ) : I afi| < 1, x 2 = -1} U {x = (x 1 ,x 2 ): \xi \ < 1, x 2 = 1}, 
We use the following boundary conditions 



(1.5) 



l! u t||i°°(0,T;L 2 ) + ||Vu|| L oo( 0T;i 2( n )) + ||V7r|| i 2( 0i T;L 2 ) 

+ l|Vwt|| L 2 (o,T;L 2 ) + \\D 2 u\\ L 2 {0T . L2) < C, 



2. PRELIMINARIES AND BASIC RESULTS 



(2.1) 




u is 2-periodic w.r.t x\. 
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Here, x\ represents the tangential direction to Y and this idealized setting of a 
"periodic strip" corresponds to the half-space, but without complications at infinity. 

Given q > 1, by L 9 (f2), we indicate the usual Lebesgue space with norm || • || g . 
Moreover, by W k ' q (n), k a non- negative integer and q as before, we denote the 
usual Sobolev space with norm || • \\k iq . We also denote by W ' 9 (O) the closure of 
Cg°(fl) in W l ' q {0) and by W' 1 ^ ' (n), q' = q/(q - 1), the dual of W 1,9 (Q) with 
norm || • Let X be a real Banach space with norm || • \\x- We will use the 

customary spaces W k ' q (0, T; X), with norm denoted by || • \\w k 'i(o,T;X)> recalling 
that W°' q {0,T;X) = L q (0,T;X). We will also use the notation fl T '■= ^ X (0,T) 
and we will not distinguish between scalar and vector fields and the symbol ( • , • ) 
will indicate a duality pairing. Here and in the sequel, we denote by C positive 
constants that may assume different values, even in the same equation. We also 
define 

V q := {v € VK M (ft) : V • v = 0, «| r = 0, v is 2-periodic w.r.t. xi}, 

with dual space V'. Since the extra-stress tensor 5* is a function not of the gradient, 
but of the deformation tensor (in order to have frame invariant equations) we recall 
a Korn-type inequality, see [5] 

Lemma 2.1. There exists a positive constant C = C(f2) such that 

IMI fl + l|Vu|| 9 <C||2>«|| 9 , for each vtV q . 

We write / ~ g, if there exist Co, C\ > such that cof < g < c\f . When 
considering the tensor S(D) = (5 + \D sym \) p - 2 D sym , introduced in (TO]) (where D 
is a second order tensor and D sym its symmetric part) it can be easily checked that 
for any second order tensor C, the following relations are verified 

2 

(2.2a) Yl d kl S l3 {D)C l3 C kl > (p - 1)(<5 + |^ m |)^ 2 |C| 2 , 

i,j,k,l— 1 

(2.2b) \d kl S^(D)\ < (3-p)(S+\D s y m \)P~ 2 . 

The symbol duSij represents the partial derivative dSij/dDu of the (i, j)-component 
of S with respect to the (k, /)-component of the underlying space of 2 x 2 matri- 
ces. Monotonicity and growth properties of S are characterized in the following 
standard lemma. 

Lemma 2.2. Assume that p £ (l,oo) and S € [0, oo). Then, for all A, B e M 2x2 
there holds 

(S{A) ~ S(B)) ■ {A sym ~ B sym ) ~ (S + \B sym \ + \A sym \)P- 2 \A sym - B sym \ 2 , 

\S{A) - S(B)\ ~ (5+ \B sym \ + \A sym \)P- 2 \A sym - B sym \, 

where the constants cq, c\ > depend only on p, and are independent of 5 > 0. 

From the elementary inequality a p < a 2 !}? -2 + 6 P , valid for all < a, < 6, and 
p € [1,2], we get the relation 

(2.4) 5% +t% ~ (5 + t)^t + 6%, S,t>0 

with constants depending only on p (see Corollary 2.19]). 

Since in the Dirichlet case we need to handle in a different way tangential and 
normal derivatives, we denote by D 2 u the set of all the second-order partial deriva- 
tives of u. In addition, the symbol D 2 u denotes all partial derivatives d 2 k Uj, except 
for the derivative d 2 2 ui, namely 

2 

\Dlu\ 2 :=\d2 2 u 2 \ 2 + Yl l^l 2 - 

i, j, k — 1 
(i,K) * (2,2) 
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We introduce the following quantities strictly related to the stress tensor S and 
coming naturally in the problem, when using the techniques introduced in pi fTTlLT^] : 

(2.5a) 2i(«):= [ (5 + iVulf^^Vufdx, 

Jn 

(2.5b) I[u):= f {5 + \Vu\y- 2 \VVu\ 2 dx, 

Jn 

(2.5c) J(u):= [ (5+\Vu\) p ' 2 \Vu t \ 2 dx, 

Jn 

where I is obtained by integration by parts when testing the extra stress-tensor S 
with — Au (and this is possible in the periodic-case); a multiple of I\ is obtained 
testing with — dull and the calculations are possible in the flat domain; Finally 
a multiple of J is obtained testing with uu and calculations are valid also in the 
Dirichlet case, for a generic domain. 

We will also use this classical result, see Necas [20] . 

Lemma 2.3. If it holds V.g = divG, for some G G (L q (n)) 2x2 , for 1 < q < +oo 
then 

"g - I g{x)dx <c\\G\\ 



i 

Let us recall the definition of weak solution to the Problem (|l.ip - (|2.1j) . 

Definition 2.1. Let T > and assume that f G L 2 (0, T; V 2 ')- We say that u is a 
weak solution of problem (|f .![) if: 

(2.6a) u G L 2 (0, T; V 2 ) n L°°(0, T; L 2 (f7)), 

(2.6b) Mt eL 2 (0,T;V 2 '), 

u(t)ipdx + vo / / V«(s) dsds + / (S(T>u(s)), T>ip} ds 
(2 6c) Jn 

/ (u(s) • V) <p u(s) dxds = / UQtpdx+ / (f(s),<p)ds V^GVa- 

Due to the fact that > 0, the existence of weak solutions follows for all p > 1 
in a standard way, and one has not to resort to very sophisticated tools as in 
Dicning, Ruzicka, and Wolf [12] ■ We will come back later on, for the motivation on 
this assumption on vq. In particular, we do not have any further restriction on p 
and the proof follows the same lines of the classical work on monotone operators, 
as summarized in Lions [18] . The result below is part of the folklore associated 
with non-Newtonian fluids. We will give a sketch of the proof since some of the 
calculations will be used many times in the sequel. 

Theorem 2.1. Let be given Uo,Vi > 0, p G [1,2], uo G L 2 (fl) with divuo = 
and (uq ■ n)|r = 0, and f G L 2 (0,T, V 2 '). Then, there exists a unique solution u 
to (|1.1[) - (|2.1[) satisfying ([2.6a|) - (|2.6cp . Moreover, the following estimates are verified 

IMIl~(o,t ; l2) + ^o|| Vm||| 2(0 ^. jL2) < C 

ll U t|li 2 (0,T;V 2 ') ^ 

where C = C{p, S,v> 0) ui, \\f\\ L 2 (o,T;V^, W^oh, T, O). 

Proof. We deduce the a priori estimates on which the existence of weak solutions 
to (|1.1|) - (]2.1|) is based. More properly, one should consider approximate Galerkin 
solutions defined as follows. Let {w r }, with r G N, be the eigenfunctions of the 
Stokes operator and let {A r } be the corresponding eigenvalues; we define X m :— 
spanjej 1 , . . . ,uj m } and P m is the orthogonal projection operator over X m . We will 
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seek approximate functions u m (t,x) = Y^T=i c™(t)uj r (x) as solutions of the system 
of equations below, for all 1 < r < to, t <E [0, T] 

u™u r + j/ Vw m Vw' r + v x S{Vu m ) Vu r + (u m • V)u m uj r ] dx = (/, uj r ), 
u m (0)=P m u o . 

Taking the L 2 -product of (ll.laj) with u m , using suitable integrations by parts and 
Young inequality we get 

+ "ol|V»'"||l + ^(i+|xv"ir J |i>»"T dx 

Using (12.41) and integrating in time we arrive at the following inequality 



\\u m (i)\\l + vo 



f \\Vu m {s)\\ 2 2 ds + C Vl f \\Vu m (s)\\Pds 
Jo Jo 



1 



< W1 + - / II/(S)II-1,2 ^ + C(»^l^, 

vo Jo 

for a.e. t e [0,T]. We estimate, by comparison, the time derivative. The only term 
which requires some care is the extra-stress tensor S. Since p < 2 we get 

(S (Vu m ),V<p) ds < 1 1 S{Vu m ) IU2 ( n T ) || V<^ || i2(nr) 

< 11^11^(^)11 v^|U w 

<C(T, 0)1^^11^ ||V^|U 2( n T) . 

Whence, by standard calculations 

(2.7) [ \\uT(s)\W 2 ds<C, 

Jo 

for a constant C depending on p, i/q, v\, ||/||l 2 (o.T;V)' ll M oli2, T, and fi. This proves 
that if u m is a Galerkin approximate solution then, uniformly in to G N, 

u m E L oo {Q,T;L 2 {n))nL 2 (0,T;V 2 ) and £ L 2 (0, T; V 2 r ). 

Note that we can extract sub-sequences converging weakly to some u in L 2 (0, T; V2), 
weakly* in L°°(0, T; L 2 (f2)) and, by Aubin-Lions theorem, strongly in L 2 (£It), and 
a.e. in Qt- We have enough regularity to pass to the limit in the convective term. 
Moreover, since S(Du m ) is bounded uniformly in L 2 (1V), it follows that S(T>u m ) — >• 
A for some A in L 2 (VLt). (Observe that without the Laplacian term we would have 
only a bound in LP (^t))- We have now to check that A — S{T>u). This is obtained 
with the monotonicity trick, see e.g. [181 §2-5.2]. By usual Sobolev embeddings 
(since we are in two dimensions) the function 1 1— > J n (u ■ V) u u dx S L 1 (0, T) , hence 
we can write the energy equality between any couple < sq < s < T 

(2.8) l -\\ u { s )\\ 2 2 + v f \\Vu\\ 2 2 dT + Vl [ S {A,u)dT=hu(s )f 2 + [\f,u)dT. 

1 J SQ J SQ 1 J SQ 

Defining for <fi £ L 2 (0, T; V2) (a test function with the same regularity of u) 

X™ :— vi f (S(Vu m ) - S(V(f>),Vu m -V^dr + vo f \\ Vu m \\ 2 d T + \\\u m {s)\\ 2 , 
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it follows, by using that S is monotone and by semi- continuity of the norm, that 
liminf m X™ > Vq ||Vw||2 dr + i || -u-(s) || i, and also that 

limA7= f\f,u) + hu \\l-vi f \A,V<I>) dr - v x f {S{Vct>),Vu - V<j>) dr. 



2 

Hence, by using the equality (|2.8I) we get 

vi [ (A - S(D<p), Vu - T>4>) dr >0 a.e. s € [0,T]. 
Jo 

We fix <j> = u - Xtjj for ip e L 2 (0, T; V 2 ) and A > 0. Finally, letting A — > 0+ the 
thesis follows. 

It is important to point out that the weak solution above constructed is unique. 
Let us suppose that we have two solutions u± and u 2 corresponding to the same 
data. We obtain the following inequality for U := U\ — u 2 (This follows by using 
the usual interpolation inequalities as for the Navier-Stokes equations and since U 
is allowed as test function, see Constantin and Foias [7]) 

\\U(t)\\ 2 2 + iso f \\VU(s)\$da + Vi ( (S{Vu l )-S{Vu 2 ),Vu 1 -Vu 2 )ds 
Jo Jo 

<- Avui( S )||l||C/( S )||^ S . 
v o Jo 

Since S is monotone (cf. Lemma \2. 2 1) the integral involving the extra stress-tensor 
is non-negative. Using the Gronwall lemma and the energy estimate one obtains 
that U ee 0. □ 

This latter result is very relevant since it allows to conclude that all the sequence 
{it m } converges to u. Moreover, if we have other a priori estimates on u m , the extra- 
regularity is inherited by weak solutions directly. This will be used in the proof of 
Theorem 11.11 Observe also that, at moment, we do not have any information on 
the pressure, apart that there exists as a distribution, by using De Rham theorem. 



3. The space-periodic case 

In this section we are concerned with the space-periodic case, that is 57 = T 2 . 
Each considered function w will satisfy w(x + 2nei) — w(x), i = 1,2, where {e±, e 2 } 
is the canonical basis of R 2 . We also require all functions to have vanishing mean 
value, to ensure the validity of the Poincare inequality. We prove some regularity 
results and we will show why the hypothesis v$ > seems necessary in many 
arguments. We define V por (f2) as the space of vector-valued functions on ft that are 
smooth, divergence-free, and space periodic with zero mean value. For 1 < q < oo 
and k € N, set 

w d£ (°) : = {closure of V pcr (ft) in W k > q (n)}, 

endowed, with the usual norms. 

In the space periodic setting many calculations are simpler since we can use 
— Au as test function (now formally but the procedure goes through the Galerkin 
approximation). Since in the 2D space-periodic case Jq(u ■ V) uAudx = we get 

(3.1) 4l|V U || 2 + ^o||A U || 2 + v x X{u) < C||/|| 2 , 

dt 

hence, if we are able to construct such a solution (this is not trivial at all due to 
some technical issues when passing to the limit in / Q X{u m {s)) ds, for a fixed T > 0) 
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that and if vq — we obtain as higher order estimate 

/ X(u)dt= [ [ (5+ \Vu\) p - 2 \VVu\ 2 dxdt < +oo. 
Jo Jo Jt 2 

We recall the following lemma, which is an adaption of [SJ Lemma 4.4] to the two 

dimensional case. 

Lemma 3.1. Let p G (1,2], S £ (0,oo), and I g [1,2). Then, for all sufficiently 
smooth functions u with vanishing mean value over fl, the following relations hold 
true 

\\u\\ p 2e <c(l{u) + SP), 

Hence, the information on the regularity in the space variable which we can 
extract from (|3.1[) . in the case v$ = 0, could be at most 

u E W 2 ' e {T 2 ) V£ < 2, a.e. t e [0, T]. 

This is not enough to employ Thm. 11.11 and explains the introduction of the hy- 
pothesis vq > 0. 

Proof of Proposition In the light of the above observations the proof follows as 
in the 2D Navier-Stokes equations, see[10]. We test the equations by — t Au m and 
we have 

j(t\\Vu m \\ 2 ) + u t\\Au m \\ 2 + ly.tliu" 1 ) < C*i||/|| 2 + ||Vu|| 2 . 

Hence, no matter of the non-negative term coming from the extra-stress tensor, 
integrating in time over [0,T] we have that Viu m € L 2 (0, T; W 2 - 2 (T 2 )). Due to 
uniqueness of the solution the whole sequence {u m } converges to u and by lower- 
senhcontinuity of the norm we obtain that y/tu G L 2 (0, T; W 2 ' 2 (T 2 )). □ 

For the sake of completeness, we recall that in the periodic 2D case, with vq = 
it is possible to prove the following result of existence of regular solutions, see 
Kost 17 , which is an adaption of those in [S] for the 3D case. (Observe that in 
absence of the Laplacian also the existence and uniqueness of weak solutions is 
more delicate and the limit process on Galerkin solutions requires some care) . The 
following result, which is of interest by itself, is not enough for our purposes of 
studying uniqueness for solutions to (|1.4p . 

Proposition 3.1. Let be given S g [0,<5o]> for some 5o > 0, set vq = 0, v\ > 
0, and let p € (1,2]. Given T > 0, assume that f € L°°(0, T; H /1,2 (T 2 )) n 
W^'^T^CIT 2 )). Let u Q G W 2 ' 2 {T 2 ) be such that divu = and div S{Vu Q ) € 
L 2 (T 2 ). Then, there is a time < T' < T (depending on the data of the prob- 
lem) such that the system (|1.1[) , has a strong solution u on [0,T"] satisfying, for 
re (4/3,2), 

ue L q (0,T';W 2 ' r (T 2 ))DC(0,T';W 1 ' q (T 2 )), Vq < oo. 

Remark 3.1. One can obtain further regularity results for Ut and also for \7ir (the 
latter if S > 0). 

4. Space- time regularity in the Dirichlet case 

In this section we consider the time evolution problem with Dirichlet boundary 
conditions and we prove a result of regularity for smooth data. Then, we will relax 
some of the assumptions to prove the main result of the paper. We start by showing 
a first regularity result for the time derivative of the solutions to the problem 
with Dirichlet boundary conditions. We prove now some results by using as test 
functions first and second order time derivatives of the velocity. These are legal test 
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functions, since if u is divergence-free and u\r = 0, then ^£ shares the same two 
properties, for all k € N. In particular, the following result is valid in any smooth 
and bounded domain, while the hypothesis of fiat boundary will be used for the 
W 2,2 (fi)-regularity. 

Lemma 4.1. Let p G (1,2], 5 > 0, / € W 1,2 (0,T; L 2 (Q)), u G W 2 ^ 2 (n) n Va, a«rf 
Ze£ u be a weak solution of problems (|1.1[) - (|2.1I) . Then, 

f4 , H u t|lL~(0,T;L2) + II Vu IIl~(0,T;L2) + ^0||VUt||l2 (0)T . i2) 

+ n||.7(u)||ii(o J T) <C, 

where the constant C depends on p, 8, vq, v\, 2(0, T:L 2 )7 || u o||2,2; T, and Q. 

As in the previous result we only prove the a priori estimates. A complete proof 
can be obtained through a Galerkin approximation and for the reminder of this 
section we drop the superscript "to" . We also define 

M{t) := I (5 + s) p - 2 s ds > 0, for t > 0. 
Jo 

Observe that M(t) ~ (S + tf- 2 t 2 and also (6 + t)P- 2 t 2 < t p , with 1 < p < 2. This 
shows that 

(4.2) 0<M{u) := / M(\Vu\)dx <C{p)\\Vu\\ p , with 1 < p < 2. 

Jn 

Proof of Lemma \4-1\ First, we multiply (|l.la|) by u t and integrate by parts. We 
observe that taking the duality of — div S(T>u) against u tl we get 

(4.3) - (div S(Vu),u t ) = (S(Vu),Vu t ) = j t M{u). 
By suitable integrations (since &\v u t = 0) we obtain 

IK||| + Y S ||Vu||| + n^M(u) = / (fu t -{u-V)uu t )dx. 

By using Holder and Gagliardo-Nirenberg inequalities, with the boundedness of the 
kinetic energy, we get, for all e > 



(u ■ V) it itt (ix 



< ||u||4||Vu|| 2 |M 



^CHulllllVulllllVullallutllJIlVutlll 
<c e (||Vii||| + ||ut||i||Vu||l) +e||Vttt||l. 
Thus, we obtain the following differential inequality 

(4.4) WutWl+j^oWVug+^Miu)) < c £ (||V U || 2 + || Ut || 2 ||V U || 2 +||/||2)+ e ||V Ut || 2 , 
which we clearly cannot use directly, due to the lack of control for Vut ■ 

Remark 4.1. Another path will be that of using improved estimates for Vm to 
estimate the convective term, see the last section. 

We take now the time derivative of (|l.la[) . multiply by u t and integrate by parts 
(recalling that Jq{u ■ V) u t u t dx = 0) to obtain 



(4.5) ~\\u t \\ 2 2 +iy \\Vu t \\ 2 +^(d t (S(Vu)),Vu t ) < 



((ut ■ V)«w 4 + ft Ut) dx 
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By (|2.2a[) the term involving S in (|4.5p is non-negative being estimated from below 
by a multiple of J{u) > 0. Let us focus on the right-hand side of (|4.5[) . By using 
Holder and interpolation inequality, and the energy estimate we get, for each 77 > 0, 



(lit • V) u Ut dx 



(u t • V) Ut u dx 



< ||ii t || 4 ||Vu t || 2 ||u| 



4 



<c|KIIJ||v Ut |||H|l||v u ||l 

< c.jllVullillutlll + TjIlVutlll, 

hence, choosing rj > small enough we get 

(4.6) j t \\u t \\l + voWVMl + < CdlVulllHutl^ + Wftg). 

Summing up (|4.4I) - (|4.6[) and choosing e > small enough we get finally 



-77 (||«t|la + «ol|V«||a + viM(uj) + \\u t \\ 2 2 + H)||Vitt||l + v x J{u) 

<c(||v«||li|tt t ||l + ||v«ni + H/HI + H/tiiS). 



To integrate over [0, T) we need to make sense to ||itt(0, -)||2- From the assumptions 
on the data, the fact that S > 0, and u m (0) — P m uo we easily get (cf. O § 5]) that 

ikr(o)ii 2 <c(iwiii, 2 +iir(o)ii 2 ). 

Recall that we are working on the finite dimensional approximation u m and taking 
the limit m — > +00. With Gronwall lemma and by using the fact that Vu S 
L 2 (0,T;L 2 {n)), we get for a.e t £ [0,T] 

Il«t(*)ll2 + ^l|Vu(t)||l+ (|K(«)||2 + ^o||Vu t (s)||| + J(u(s))) ds 

< C(vq,Vx,6,T, \\f\\w^{0,T;L2), \\ «o|| 2,2, H) , 

hence the thesis. □ 

Remark 4.2. The hypotheses on the external force can be slightly relaxed, but this 
is inessential in our treatment. 

We now prove Proposition 11.21 For the reader's convenience we split the proof 
into two parts. First, we perform a preliminary study of the system obtained 
removing the convective term (u • V) u from (jl.lal) . 

(4.7a) u t - i/ Au - v x div S{Vu) + Vtt = / in [0, T] x Q, 

(4.7b) divu = in [0,T] x Q, 

(4.7c) u = in [0,T]xT, 

(4.7d) u(0) = M in O, 

and focusing on the role of the nonlinear stress-tensor. The system (|4.7I) can be 
treated similarly to a steady state problem if we have good enough a priori esti- 
mates on u t . We will then address the full problem fj 1 . 1 [) - (12 . 1[) . by adding suitable 
estimates for the convective term. 

Lemma 4.2. Let vq > 0, 8 > and p £ [1)2]- Given T > 0, assume that 
u e W 2 ' 2 (Q) n F 2 and / e ^^(O.TjL 2 ^)). I7ien, proWem (|477 |l - <P^|I admite a 
unique solution, such that (|1.5[) ZioWs frwe. 



Proof. We adapt to the time-dependent case a technique with three intermediate 
steps taken from [51 |5J: In the first step we bound the tangential derivative of 
velocity and pressure; In the second step we estimate the normal derivative of the 
velocity field; In the last step we estimate the normal derivative of the pressure. 



UNIQUENESS OF PARTICLE TRAJECTORIES FOR NON-NEWTONIAN 2D FLUIDS 11 



Again we merely prove the a priori estimates. Observe that for this simpler 
problem without convection, the same existence proved in Theorems 12.11 and reg- 
ularity from Lemma 14.11 clearly hold true (this is particularly relevant for what 
concerns Ut). 

Step 1. We first prove that the following estimates, concerning the tangential 
derivatives, hold true 

(4.8) ^ ||V9im|| 2 + ^o||<9 2 2 2 M 2 ||2 2(OT;i2) + ||5i7r||£ a( o,T i £») < C, 

where C depends on p, 5, u , vi, \\f\\ w 1 - 2 ^^-^ 2 ), ||"o||2,2, T, and fl 

We now use the particular features of the flat domain. Multiplying equation ()4.7a[) 
by —d\ x u and integrating by parts, it follows that 

^IMl + ^ol|va lU || 2 + ( P - i)vif v + \vu\y- 2 \d<Du\ 2 dx < \\f\\ 2 \\d 2 llU \\ 2 . 

By applying Young inequality and using relation (I2.5a|) . we get a.e. in [0,T] 
\\d lU (t)\\ 2 2 + f ( Vo \\Vd 1 u{s)\\ 2 2 + v 1 l 1 {u{s)))ds 



(4.9) 



<C(||V U0 || 2 + - f ||/( 
v a Jo 



sW 2 ds 



and, since div it = 0, d 22 u 2 = —d 21 ui and the estimate on d 22 u 2 follows. 

Let us focus on the pressure term. Differentiating the equation (|4.7ap with 
respect to the tangential direction xi, one has that 

V9i7r = Vq div d\ Vu + v\ div d\ \{5 + \T>u\) p ~ 2 T>u\ — diUt + dif a.e. in Qt- 

We observe that diu t — div \ ^ Ul l \ and dif — div (y |? ]. Hence to apply 

\o t u 2 Uy \j 2 Uy 

Lemma |2~3"1 to estimate din, we only have to bound the term di [(6 + \T>u\) p ~ 2 'Du\ . 

A direct computation gives 

T>ii 

d x [(S + \Vu\) p - 2 Vu] =(S+ \Vu\) p - 2 diVu + (p - 2)(S + \Vu\) p ' 3 {Vu ■ diVu)- 



\Vu\ 

and consequently 

\di [{5 + \Vu\) p - 2 Vu] | < (3 - p)(8 + \Vu\) p - 2 \diVu\ a.e. in Q T . 
Therefore, by comparison di [(5 + \Du\) p ~ 2 T>u\ g L 2 (£l) and it follows that 

f \di[(S+ \Vu\) p - 2 Vu}\ 2 dx < c5 p - 2 li{u) a.e.t£[0,T]. 

By applying Lemma 12.31 we have that 

\\dnr\\ 2 2 < \\u t \\ 2 2 + u \\diVu\\ 2 + viC8 p - 2 Xiiu) + \\f\\j a.e. t E [0,T], 

from which, integrating in time over [0, T], using (|4.9[) and recalling the bounds 
previously proved on u t , diVu, and Xi, then (|4.8I) follows. 

Step 2. To bound ||c?f 2 ' u i lli 2 (o,T:Z, 2 ): we consider a narrower range of values for 
the parameter p. Under the same hypotheses as before, but for p G [§,2), we have 

II<9 2 2 w iIIl 2 (o,T;L 2 ) < C, 

where the constant C depends on p, 5, v$, vi, \\fWw 1 - 2 (o, t ; l 2 ), 1 1 1 1 2,2 , T, and £1. 

We follow the main lines established in the proof of [51 Lemma 3.3]. By calcu- 
lating d 2 [(5 + \T>u\) p ~ 2 T>u] , the first equation in (|4.7a[) can be written as 

(4.10) aid\ 2 ui = -Fi - f x + d t ui + dnr, 
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where 

ai := u + ^ (S + \Vu\f- 2 +^{p- 2) {& + {Vu) 12 {Vu), 



and 



Fx := [v + vi(t + \Vu\) p - 2 ]d* lUl + ^{6 + \Vu\y- 2 dl 2 u 2 



(8 + \Vu\) 



p-3 



\Vu\ 



k,l = l 



+ vi(p-2) 

By direct calculations it can be easily seen that 



[ (Vujud^VujutVu)^ + dl 2 u 2 (Vu) 22 {Vu) 12 



d\ 2 u 2 



(Vu) 12 {Vu) 



m < c 



+ "l (p ~ ~) V + \Vu\y- 2 ] \Dlu\ a.e. in Q T 



and by using that p > | we get 



u + ^{p-^)(5 + \vu\y 



> vq > 0. 



Division of both sides of (|4.10j) by ct\ is then legitimate and we infer that 

\dl 2Ul \ < C(\D 2 u\ + — (|Si7r| + \d tUl \ + a.e. in Q T . 

Therefore, squaring and integrating over Ot we get 

||&«i(s)||!<k < - f T (\\D 2 u(s)\\l + \\dxir(s)f 2 + \\d tUl ( s )\\l + \\h(s)\\t) ds, 
Jo \ / 

which, by the previous results is finite. This finally shows that D 2 u g L 2 (Q,t). 

Step 3. The final step, which is not strictly required for the particle trajectories 
uniqueness, is the regularity of the normal derivative of pressure. Nevertheless, we 
include it for the sake of completeness. Under the same hypotheses we have 

\\d 2 ir\\L 2 (0,T;L 2 ) < Cj 

where the constant C depends on p, S, v<y, u\, \\f\\ wi. 2 (o, T;L 2 ), 1 1 i^o 1 1 2,2 , T, and il. 
By using the second equation in (I4.7a[) . one can write 

|9 2 7r| < c(vq + MP- 2)(S + \Vu\f- 2 ^ \D 2 u\ + \8 2 u t \ + \f 2 \ a.e. in Q T . 

Hence, straightforward calculations lead to 

i-T 

\\d 2 TT(s)\\ 2 d S < C 



\\D 2 u(s)\\ 2 2 + \\d 2 u t ( S )\\i + \\f( S )\\ 2 2 )d S 



and the assertion follows as a consequence of the previous results. 



□ 



We finally prove the same regularity results also in presence of the convective 
term. We use a perturbation argument, treating (it • V) u as a right-hand side in 
equation (|l.la|) . 



Proof of Proposition Here, we use the a priori estimates obtained for the prob- 
lem (|4.7I) with external body force 

T ':= -(u- V)it + /. 

In the derivation of estimates for it t we used that || /||w/i, 2 (o, t ; l 2 ) > while in Lemma l4.2 
the estimates depend essentially on the L 2 (r2r)-norm of the external force. Hence, 
by using Lemma 14.11 it is then sufficient to estimate \\(u ■ V) u|jx, 2 (o,T;i 2 ) m terms 
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of second order derivatives of u, to follow the same calculations in Step 1-3 of the 
previous result. 

By applying Holder, Gagliardo-Nirenberg, and Young inequalities and the energy 
estimate, we get for each e > 



By using the same calculations as in the previous proposition and the a-priori 
estimates in (|4.ip -especially that V« € L°°(0, T; L 2 (f2))- we have 



As a consequence of the above result we have full L 2 -space-time regularity of 
the solution up to second-order space-derivatives, hence the uniqueness of particle 
trajectories. The result is not optimal in view of application to uniqueness of 
trajectories, in the sense that some of the hypotheses can be slightly relaxed. For 
instance f t £ L 2 (Qt) and u e W 2 ' 2 (Q) can be removed (at the price of less 
regularity on u t ) by following a slightly different path as we do in the next section. 



In this section we finally address the problem of the uniqueness of particle tra- 
jectories under "minimal" assumptions on the data. We will show how the previous 
regularity result, together with Theorem 1 allow us to prove Theorem 11.21 

Proof of Theorem In the same way as in the proof of Lemma 14.21 we perform 
separately the a priori estimates for the normal and tangential derivative of the 
time- weighted \fl u m (which we call y/tu). In particular, here we do not use a lot of 
regularity on ut, but we have to face with a non-smooth uq. By adapting standard 
weighted estimates, we multiply the equation (|l.la|l by —td^u. Integrating by 
parts, and with Young inequality we obtain 



i i 



(4.11) 



||(u-V)u|| 2 <||n||4||Vu||4<c|| U ||I||Vu||I||V U ||I||£' 2 u| 

< c £ \\Wu\\l + e\\D 2 u\\ 2 . 




and, by choosing e > small enough, we end the proof. 



□ 



5. Proof of Theorem 11.21 




Integrating in time and using the energy estimate to bound J Q ||<9iu|| 



2 ds it follows 



(5.1) 
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We take now the L 2 -inner product of (|l.la[) with t ut . By suitable integrations by 
parts, and using (|4.2[) - (|4.3[) we reach 



iWi + y|(t||V«||l)+^(iM(«)) 



dt 



< t 



(u ■ V)ttut dx 



f u t dx 



+ iy \\Vu\\ 2 2 + Cis 1 M(u) 



t 



< i (ll(" • V) u\\ 2 + \\f\\l) + -\\u t \\ 2 + u \\Vuf 2 + C^\\Vu\\p. 

Integrating this inequality in time, by appealing to the energy inequality and re- 
calling that M(u) > 0, it follows that for a.e. t e [0,T] 



u t || Vt*(t)||| + / s \\u t (s)\\l ds 



(5.2) 



< C 



\\u \\l + / s(\\(u(s)-V)u(s)\\l + \\f( S )g)ds 



where C depends on p, S, vq, v\, T, and SI. 

Let us now focus on the normal derivatives of u. Arguing as in Step 2 of the 
proof of Lemma H~2l and replacing / with / + (u ■ V)u, we infer that 



\d 2 22Ul \ <c(\Dl 



2v 



|Qi7r| + |«t| + |(u-V)u| + |/| 



a.e. in fix- 



Then, squaring, multiplying by t, and integrating over (0, t) x SI, we find 



(5.3) 



s\\dl 2Ul {s)f ds 



< 



C 



v o Jo 



(\\Diu(s)\\t + \\dMs)\\l + \\ut(s)\\l + \\(u(s) ■ V) u(s)\\l) ds, 



To control J* s ||9i7r(s)||| ds we proceed again as in Step 2 of the proof of Lemma l4.2l 
Thus, for a.e. t € [0,T], the following inequality holds true 



s 1 1 c?i7r(s) 1 1 1 ds 



. (\\u t (s)\\l + ||axVu(*)||i + SV-^uKs) + \\f(s)\\l + \\(u(s) ■ V) u(s)\\ 2 ) ds 
u \\l+ f s(\\{u{s)-V)u{s)\\l + \\f{s)\\l)ds , 



< C 

< c 



where we have used relations (|5.1[) and (|5.2I) . Once again we apply (|5.1|) . so that 
relation (|5.3I) gives, for a.e. t £ [0, T] 



s\\dl 2Ul (s)\\ 2 2 ds<C \\u Q f 2 + / s(\\(u(s)-V)u(s)\\i + \\f(s)\\ 2 2 )ds 



where C depends on p, 6, i/q, v\, T, and Q. Summing up the above inequality 
with (jSTTj) and (pT2"j) . we get for a.e. t G [0,T] 

t||Vu(t)||| + i/ / s\\D 2 u(s)\\ 2 2 ds 



< C 



KHI+ / s(\\(u{s)-V)u(s)\\ 2 2 + \\f(s)\\l)ds 
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whit C depending on p, 6, i/q, v\, T, and Q. The convective term can be estimated 
as in (|4.1ip and choosing e > small enough we get, for a.e. t £ [0, T], 

t\\Vu(t)\\l+v [ s\\D 2 u(s)\\lds 
Jo 

<c £ f (s HVu(a)Hi) || Vu(s)f 2 ds + C{ Pl 6, uo, u 1} \\f\\ L ^ T . m , ||uo|| a , T, (I). 
Jo 

Hence, by using Gronwall inequality over [A, T] (for any A > 0), letting A — > + , 
and by using the energy inequality we get 

/ t\\D 2 u(t)gdt< C(j>,8,V ,V U ||/|| i a ( 0,T i £»),||«0||2,r > n). 

Jo 

Then, the assertion follows by means of Theorem 11.11 □ 
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